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Abstract. The one dimensional symmetric simple exclusion process (SSEP) is one of 
the very few exactly soluble models of non-equilibrium statistical physics. It describes 
a system of particles which diffuse with hard core repulsion on a one dimensional 
lattice in contact with two reservoirs of particles at unequal densities. The goal of this 
note is to review the two main approaches which lead to the exact expression of the 
large deviation functional of the density of the SSEP in its steady state: a microscopic 
approach (based on the matrix product ansatz and an additivity property) and a 
macroscopic approach (based on the macroscopic fluctuation theory of Bertini, De 
Sole, Gabrielli, Jona-Lasinio and Landim). 

PACS numbers: 02.50.-r, 05.40.-a, 05.70 Ln, 82-20-w 
1. The Symmetric Simple Exclusion Process 

Understanding the steady state properties of systems in contact with two heat baths 
at unequal temperatures or two reservoirs of particles at unequal densities is a central 
question in the theory of non-equilibrium systems [U El EJ H]. Here I would like to 
focus on one exact result which was obtained during the last decade on the steady 
state of one of the simplest models of a non-equilibrium system, the one dimensional 
symmetric simple exclusion process (SSEP). For the SSEP this exact result, which gives 
an expression of the large deviation functional of the density, can be derived either from 
a microscopic description of the steady state [HI El [7j or from a macroscopic approach, 
which was developed by Bertini, De Sole, Gabrielli, Jona-Lasinio and Landim [8| [9| [TO]. 
These two approaches are discussed below. 

The symmetric simple exclusion process (SSEP) describes a gas of particles diffusing 
on a lattice with an exclusion rule which prevents a particle to move to a site already 
occupied by another particle [HI [121 [IS [12 US] . Here we consider the one dimensional 
version with open boundaries. The lattice consists of L sites, each site being either 
occupied by a single particle or empty. During every infinitesimal time interval dt, each 
particle has a probability dt of jumping to the left if the neighboring site on its left is 
empty, dt of jumping to the right if the neighboring site on its right is empty. At the two 




Figure 1. The symmetric simple exclusion process: particles diffuse with hard core 
repulsion on a one dimensional lattice connected at its ends to two reservoirs of particles 
at densities pa and pi,- 



boundaries the dynamics is modified to mimic the couphng with reservoirs of particles 
at densities pa for the left reservoir and pb for the right reservoir: at the left boundary, 
during each time interval dt, a particle is injected on site 1 with probability 2padt (if 
this site is empty) and a particle is removed from site 1 with probability 2(1 — pa)dt (if 
this site is occupied). Similarly on site L, particles are injected at rate 2pft and removed 
at rate 2(1-/9;,). (The factors 2 in the boundary rates simplify some expressions below 
but do not affect the large scale properties). 

The SSEP is obviously a model of transport of particles between two reservoirs at 
densities pa and ph- It is also a simple model of heat transport between two heat baths 
at temperatures and T^, if one interprets each particle as a quantum of energy e, with 
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Under the evolution rules of the SSEP, the system reaches, in the long time limit, a 
steady state. If one divides the system of length L into n boxes of size / (with of course 
L = nl), one can try to determine, in this steady state, the probability of a certain 
density profile {ri,r2..r„}, i.e. the probability of seeing Iri particles in the first box, lr2 
particles in the second box, ... /r„ particles in the nth box. For large L, one expects the 
following L dependence of this probability 

ProL(ri, ...rn\pa, Pb) ~ exp[-LJ'„(ri, r2, ...r„|pa, pb)] (2) 

where J-'n(yi^i,i^2, ■■■^n\Pa, Pb) is called the large deviation function [TB] of the density 
profile {ri,r2..r„}. When the number n of boxes becomes large, keeping the number / 
of sites in each box also large, one can introduce a continuous variable x = k/n, the 
densities ri,r2..r„ become a density profile 
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p(x) = p{- 
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and the large deviation function v^(ri, r2, ...r„|pa, Pb) becomes a functional 
J^{{p{x)}\pa, Ph) of the density profile p{x) 

Vio l{{p{x)]\p a, pb) ~ ew[-LF{{p{x)]\pa, Pb)] ■ (3) 
The result which was obtained in [5l [9l |6] is that, in the non-equilibrium case, i.e. for 
Pa 7^ Pfe, the exact expression of the large deviation function J^{{p{x)}\pa, Ph) is given 
by 



•^({p(^)}|Pa,Pb) 
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where B{p, r) is given by 
5(p,r) = (l 



P) lo^ 



1-p 



plog 



p 



and where the function F{x) is the monotone solution of the differential equation 
, , ^ F(1-F)F" 

satisfying the boundary conditions -F(O) = pa and F{1) = pb. 
2. Properties of this large deviation function 
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Before describing the two main approaches which led to ( Ipp ), let us discuss briefly a 
few properties of the functional J-". 

One can solve perturbatively ^ for pa close to pb and get 

{pa - PbY 



Pa - (Pa - Pb)x - 
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y (piy) - Pa) dy 



pail- pa) 

+x (1 - y) {p{y) - Pa) dy 
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Therefore in the limit pb Pa, the expression ^ reduces to 
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+ O [{pa - Pb)' 



J^i{pix)}\pa,Pa) 



dx [B{p{x),pa)] 
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This is not surprising as, when pa = Pb, the dynamics satisfies detailed balance and 
in the steady state all the lattice sites are occupied independently with probability pa 
(Bernoulli measure). In this equilibrium case J-" is a local functional ([s]) of the density 
profile p{x). This is a special case of the much more general fact [TJ [I7j that, for any 
system (with short range interactions) at equilibrium, in contact with one or several 

reservoirs at density pa, the functional J-" is always local and is given by 

»i 



J^{{p{x)}\pa,pa) 



dx 



fiPi^)) - fiPa) - (P - Pa)f{pa) 



(9) 



where /(p) is the free energy per unit volume at density p, defined as /(p) = 
— \\mv^oo^ogZ{y,Vp)/V where Z{y,N) is the partition function of a system of 
particles at equilibrium in a volume V . 

In the Pa 7^ Pb case, the large deviation functional F can therefore be thought as a 
possible generalization of the concept of free energy to non-equilibrium systems. 

As soon as the system is out of equilibrium (pf, ^ pa) the large deviation functional 
J-" becomes non-local. This is already visible in the expansion of J-" in powers of pa — pb, 
obtained by replacing F{x) by its expansion (0 into (§ 



^{{P{x)}\pa, pb) 



dxB{p{x), p* (x)) 



I (Pa - Pb)^ 
[Pa(l - Pa)Y 

where the average profile p*{x) is given by 

P*{x) = (1 - X)pa + Xpb- 



dx / dy x(l - y){p{x) - p*{x)){p{y) - p*{y)) 



(10) 
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The knowledge of the functional T allows one to determine all the the correlation 
functions: if one defines the generating function Q of the density by 



exp [L^({a(a;)}|pa,Pfe)] = (exp 



L \ a{x)p(x)dx 



(12) 



where a{x) is an arbitrary function and (.) denotes an average over the profile p{x) in 
the steady state, it is clear from (|3| that, for large L, Q is the Legendre transform of J-' 



g{{a{x)}\pa,pb) = max 



a{x)p{x)dx - J^{{p{x)}\pa,Pb) 



(13) 



By taking derivatives of (12) with respect to a{x) one can then get all the correlation 
functions. In particular 

1 s'g 



{p{x)p{y))c = {p{x)p{y)) - {p{x)){p{y)) 



L 6a{x) Sa{y) 



(14) 



a{x)=0 



A direct consequence of ( 12 ) is that all the fc-point connected correlation functions are 
long ranged and scale like L^~'^ (see [IHl [191 E!)- 



3. The steady state 

From the definition of the SSEP, if is a binary variable with Xj = 1 when site i is 
occupied and = when it is empty, one can write the time evolution of the average 
occupation (tj) 



dt 

d{n) 

dt 
dt 



2pa - 3(ri) + (rs) 
(r,_i) - 2(7-,) + (r,+i) for 2 < z < L - 1 (15) 
{rL-i)-3{TL)+2pb . 

d{n) 



The steady state density profile (obtained by writing that = 0) is [6] 

{n) = ^ (16) 



For large L, with i = Lx, one recovers the average density profile (11). 

In a similar way one can then write down the equations which govern the time 
evolution of the two point function or higher correlations. 



4. The matrix ansatz for the SSEP 

For the SSEP, one can then write down the steady state equations satisfied by higher and 
higher correlation functions, but solving these equations becomes quickly complicated. 

The matrix ansatz [201 EH [22] gives an algebraic way of calculating exactly the 
weights of all the configurations in the steady state: in [20] it was shown that, in the 
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Pro({ri,r2, ...tl}) 



(17) 



steady state, the probability of a microscopic configuration {ri, T2, ...tl} can be written 
as tlie matrix element of a product of L matrices 

{pa\X,X2...XL\pb) 

{pa\{D + E)^\p,) 

where the matrix Xj depends on the occupation number 

= nD + (1 - n)E , (18) 

and the matrices D, E and the vectors {pa\, \pb) satisfy the following algebraic rules 
DE-ED = D + E 

{Pa\2[p,E-{l-p,)D] = {pa\ (19) 

2 [{1 - p,)D - p,E] W) = \pt). 
A priori one should construct the matrices D and E (which might be infinite- 



dimensional) and the vectors {pa\ and \pb) satisfying (19) to calculate the weights of 
the microscopic configurations. However these weights do not depend on the particular 
representation chosen and can be calculated directly [201 E] from (19). 



One can calculate, using (19), the average density profile 

{pa\{D + Ey~'D{D + E)'^-^\pb) 



{p^{D + E)L\pb) 



as well as all the correlation functions and recover (16). One can also show that 

Ef\Pb) 



{Pa\{D 



L\ 



(20) 



(21) 



{pa\pb) {pa- Ph)^ ' 

(This formula is easy to derive by noticing that, for a system of size L, the average 
current is given, according to (16), by (tj — Tj+i) = (p^ — Pb)/L but is also given. 



according to ( [20|T9| ), by the ratio {pa\iD + Ef-^\pb)/{pa\{D + E)^\pb)). 



5. Additivity 



As in (17) the weight of each configuration is written as the matrix element of a product 



of L matrices, one can try to insert at a position Li a complete basis in order to relate 
the properties of a lattice of L sites to those of two subsystems of sizes Li and L — Li. 
If one defines, for arbitrary p, left and right vectors (p| and |p), which satisfy 

(p| 2 [pE - (1 - p)D] = (p| 

2[{l~p)D-pE] \p) = \p) (22) 
(note that in general (p|p') 7^ 0), it is possible to show, using DE — ED = D + E as in 



(|19|) and the property (|22|), that for pb < Pa 

dp 



{Pa\YlY2\pb) 
{Pa\pb) 



(Pa-Pb) {palYllp) {p\y2\Pb) 



Pb<\p\<Pa 



2m {pa - p){p - pb) {pa\p) {p\pb) 



(23) 



where Yi and Y2 are arbitrary polynomials of matrices D and E. (To prove (23), one 
can first prove it, using (21 ), for Yi of the form [paE — (1 — Pa)D]'^^ [D + E]'^'^ and Y2 of 
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the form [D + Ef^ [(1 - pi,)D - pi^Ef''. Then one can show, using DE - ED = D + E, 
that any polynomial Yi or Y2 can be reduced to a finite sum of such terms). 
By choosing for Yi the sum over the weights of all configurations with Iri occupied sites 
in the first a box, ... /r^ occupied sites in the kth box, and for Y2 the sum over all 
configurations with Ir^+i occupied sites in the k + 1th box, ... /r„ occupied sites in the 
nth box, one can show, using ( [l7|21[23[ ) that 

dp Li! (L-Li)! (Pa-Pb)^+' 
2z7r 



ProL(ri, ...rn\pa,Pb) 



(24) 



Pb<\p\<Pa 



Pro L^{r I,... rk\pa,p) x PrOi^L, (r^+i, ...r„|p, p?,) 

where Li = kl. This formula, which is exact for arbitrary system sizes, relates the 
properties of two disconnected subsystems of sizes Li and L — Li to those of a single 
system of size L. 



If Li = Lx, one then gets ^ for large L 

xJ'kiri, ---rklPa, F) + {1- x)J'n-k{rk+i, •••r„|F, pb) 



J^niri,r2,...rn\pa,pb)= max 

pt<F<p. 



+ xlog 



X 



x) log 



1 



X 



- \og{pa - pb) (25) 



which follows from (24) by a saddle point method (as in (24) the integration contour is 
perpendicular to the real axis, the value F of p which maximizes the integrand along 
the contour becomes a minimum as p varies along the real axis). If one repeats the same 
procedure n times, one gets 

J^n{ri,r2,...rn\pa,Pb)= ^ max -f^ J^i{rk\Fk-i, Fk)+\og ( ^'^^ ) (26) 

Pa=Fo>..>F^:>..>F„=pb ^ \ Pa — Pb J 

For large n, as F^ is monotone, the difference Fk-i — Fk has to be small for almost all k 
and one can replace J-'i(rfc|Ffc_i, F^) by its equilibrium value J-'i(rfc|Ffc, F^) = B{rk,Fk) 
(see (|8])). If one the writes F^ as a function oi k/n 



n 



(26) becomes 



^{{p{^)]\pa,Pb) =max / dx 

F(x) Jo 



B{p{x),F{x)) +\og 



F'{x) ' 

Pb - Pa 



(27) 



(28) 



where the maximun is over all the monotone functions F{x) which satisfy F{0) = pa 
and F{1) = pb- Writing the equation satified by the optimal F{x) in (|28|) leads to (g 
and this completes the derivation of ( 4p ). 
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6. The macroscopic fluctuation theory 



For a general diffusive one dimensional system of length L, in contact with a left reservoir 
at density pa and a right reservoir at density pb, the average current and the fluctuations 
of this current near equilibrium can be characterized by two quantities -D(p) and a(p) 
defined by 



lim 



lim 



m Dip) 



t 

m 



L 



{Pa - ph) for {pa - Ph) small 



for Pa = pb 



(29) 



(30) 



t L 

where Qt is the total number of particles transferred from the left reservoir to the system 
during time t. Starting from the hydrodynamic large deviation theory [231 [ISl E] 
Bertini, De Sole, Gabrielli, Jona-Lasinio and Landim [H |9l |T0] have developed a general 
approach, the macroscopic fluctuation theory, to calculate the large deviation functional 
J-' of the density ^ in the non-equilibrium steady state of a diffusive system in contact 
with two reservoirs as in figure [1} 

For diffusive systems (such as the SSEP), the density pi(t) near site i at time t and 
the total fiux Qi{t) fiowing through position i between time and time t are, for a large 
system of size L and for times of order L^, scaling functions of the form 

Pi(t)=p[i,^) , and gi(t) = LQ (31) 

(Note that, due to the conservation of the number of particles pi{t) — pi{0) = Qiit) — 
Qi+i(t), the scaling form of pi(t) implies the scaling form of Qiit)). If one introduces 
the instantaneous (rescaled) current defined by 

dQ{x,T) 



dr 



(32) 



(33) 



the conservation of the number of particles implies that 

dp{x,T) d'^Q{x,T) dj{x,T) 
dr drdx dx 

The macroscopic fluctuation theory [SI El [ID] starts from the probability of observing 
a certain density profile p{x, t) and current profile j (x, r) over the rescaled time interval 

Ti < T <T2 

2- 



;({p(a;,r),j(a;,T)}) 



exp 



-L 



T2 



dr' / dx 



jix,T')+Dipix,T')) 



dp{x,T') 
dx 



2a(p(a;,r')) 



(34) 



where the current j{x, s) is related to the density profile p{x, s) by the conservation law 
(33) and the functions D{p) and cr(p) are defined by ( 29p0 ). The physical meaning of 



(34) is that the system is locally close to equilibrium and that the fiuctuations of the 



local currents are Gaussian with averages and variances given by ( |29p0| . 

Then to calculate the probability of observing a density profile p(x) in the steady 
state, at time r, one has to find how this profile is produced. For large L, this probability 
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(|34| is dominated by the optimal path {p{x, s),j{x, s)} for — oo < s < r in the space of 
density and current profiles and 



FTOL{{p{x)}\pa, pb) ^ _ max P-oo,r ({p(x, s), s)}) 

{p{x,s),j(x,s)} 



(35) 



which goes from the average steady state profile p*{x) (given by (11) for the SSEP) to 
the desired profile p{x) 

p{x,—oo) = p*{x) ; p{x,t) = p{x) . (36) 

This means that the functional J-' of the density ^ is given by 

'j{x,T')+Dipix,T')y 



H{p{x)}\pa, Pb) 



mm 



{p{x,s),j{x,s)}J-oo 



dr' / dx 



^i\\ dp{x,T') 
dx 



2a{p{x,T')) 



(37) 



where the density and the current profiles satisfy the conservation law (33) and the 
boundary conditions (36). 



Finding this optimal path p{x,s),j{x,s) with the boundary conditions (36) is 
usually a hard problem. Bertini et al [Hj were however able to write an equation 
satisfied by J-": as (37) does not depend on r (because the probability of producing 
a certain deviation p{x) in the steady state does not depend on the time r at which this 
deviation occurs), one can isolate in the integral (37) the contribution of the last time 

(38) 



interval (r — 6t, t) and (37) becomes 



^({P(^)}) 



mm 

Sp{x),j{x) 



T{{p{x) - 6p{x)}) + 6t f dx 

Jo 



[j{x) + D{p{x))p'{x)y 



2a{p{x)) 

where p{x) — 6p{x) = p{x, r — dr) and j{x) = j{x, r). Then if one defines U{x) by 

5^({p(x)}) 



U{x) 



6p{x) 



(39) 



^^dr, one should have according to (38) 



and one uses the conservation law 6p{x) = 
that the optimal current j{x) is given by 

j{x) = -D{p{x))p'{x) + a{p{x))U\x) . (40) 

Therefore "starting" with p{x,t) = p{x) at time r and using the time evolution (for 
—oo < s < t) 



dp{x, s) 
ds 



dj{x,s) 
dx 



(41) 



with j related to p as in (40) one should get the whole time dependent optimal profile 
p{x, s) which converges to p*{x) in the limit s — > — oo. The problem of course is that J-" 
is in general not known and so is U{x) defined in (39). 

One can write from (38) (after an integration by parts and using the fact that 
U{0) = f/(l) = if p(0) = Pa and p(l) = pb) the equation satisfied by U'{x) 



dx 



'Dp' 



a 



U' 



'Dp'' 



a 



a 
2 







(42) 
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which is the Hamilton- Jacobi equation of Bertini et al [8]. For general D{p) and cr(p) 
one does not know how to find the solution U'{x) of (42) for an arbitrary p(x) and thus 
one does not know how to get a more explicit expression of the large deviation function 
^({p(x)}). 

One can however check rather easily whether a given expression of T satisfies (42) 
since U'{x) can be calculated from ( p9) ). For the SSEP one gets from ( 28p9 ) 



U{x) = log 



(43) 



:i-p(x))F(x) 

with F{x) related to p{x) by (|6|. One can then check that (42) is indeed satisfied using 
the known expressions of D = 1 and a = 2p(l — p) for the SSEP [7J (using the fact 
that F"(0) = F"{1) = which is a consequence of the fact that p(0) = F{0) = pa, 
p(l) = F(l) = pfc, and of (§). 

In fact when J-" is known, one can obtain the whole optimal path p(x, s) from the 
evolution (41) with j related to p by (40) which becomes for the SSEP 

dp{x, s 



3[x,s} 



dx 



+ a(p(x, s)) log 



p{x,s){l - F{x,s)) 
(1 -p(x,s))F(x,s) 



(44) 



where F is related to p by For (4]|6) to coincide with (37), the optimal profile p 
evolving according to (41) should converge to p*(x) as s — )■ — oo. One can check that 
this evolution (41) of p(x, s) for this current (44) is equivalent to the following evolution 
[9] of F 

dF{x, s) _ d'^F{x, s) 



ds 



dx'^ 



(45) 



where F is related to p by (|6j). Clearly (45) is a diffusion equation. Because F(0) = pa, 
F(l) = ph and because of the minus sign in (45), F{x, s) — > p*(x) as s — )• — oo. Therefore, 
due to ([6]), the density p{x, s) — )■ p*{x) as s — t- — oo. Thus ( 41p4 ) do give the optimal 
path in (37) with the right boundary conditions (36) and (37) coincides for the SSEP 
with the prediction (4][6) of the matrix approach. From ( 45||6 ) one can show that the 
time evolution of a deviation p(x, s), when it is produced is given, for small pa — Pb, by 



dp{x, s) d'^p{x, s) {pa - pb) 



'l-2p{x,s)) 



dp{x, s) 



+ O {{pa 



pb) 



(46) 



ds dx"^ Pa(l — Pa) ^ '"^ ' ds 

One can notice that as soon as pa 7^ Pb this is not the time reversal of the way a deviation 
relaxes (15) 



dp{x,t) d^p{x,t) 



(47) 



dt dx"^ 

This again is not a surprise as for non-equilibrium systems (pa 7^ Pb)i the way a deviation 
is produced (46) has no reason to be the time reversal of the way it relaxes (47). 
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7. Conclusion 



In addition to the two approaches discussed above to obtain (|4j|6]), Tailleur Kurchan and 
Lecomte [21] have developed a third approach based on a non-local change of variables 
which allows them to map the dynamics of the non-equilibrium case (p^ 7^ pb) onto the 
dynamics of the equilibrium case {pa = Pb)- 

Apart from the SSEP (and zero range processes for which the steady state measure 
is a product measure), the large deviation function J-" has been determined so far only for 
few other cases: the Kipnis Marchioro Presutti model [271 [28] . the weakly asymmetric 
exlcusion process [23 [2S], the ABC model [211 ED] on a ring for equal densities of the 
three species, driven systems [31] in particular the asymmetric exclusion process [321 ES] . 

An open question is whether one could use the macroscopic fluctuation theory to 
find the large deviation functional J-" for more general diffusive systems characterized 
by arbitrary functions D{p) and a{p) defined in ( 29p0 ). 

More recently, the macroscopic fluctuation theory has become a very powerful tool 
to calculate the large deviation function of the current in the non-equilibrium steady 
state of diffusive systems [SI ESI ESI EH EHl EHl SOI iH 112] • On the other hand exact 
calculations of the current fluctuations, starting from a microscopic model, are still very 
difficult to do [m HH Hni US] . What the large deviation functional of the density looks 
like, for a diffusive system, when conditioned on the current, remains an open question 

m- 

Looking, by a macroscopic or a microscopic approach, at diffusive systems with an 
initial condition which is not a steady state as in [HI 09], would be another interesting 
direction to pursue. 

Lastly, one knows [21 [50] that mechanical systems which conserve momentum 
exhibit an anomalous Fourier's law in one dimension. What the large deviation functions 
of the current or of the density become for such systems looks to me another interesting 
and challenging question. 
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